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\^ • We consider the anisotropic quantum Heisenberg antiferromagnet (with 

in : 

^~~^ • anisotropy A) on a square lattice using a Chern-Simons (or Wigner-Jordan) 

_!. ' approach. We show that the Average Field Approximation (AFA) yields a 

"^ ' phase diagram with two phases: a Neel state for A > Ac and a flux phase 

Q^' for A < Ac separated by a second order transition at Ac < 1. We show that 

1-^ ' this phase diagram does not describe the XY regime of the antiferromag- 

r'"! ' net. Fluctuations around the AFA induce relevant operators which yield the 

j^ ■ correct phase diagram. We find an equivalence between the antiferromagnet 

and a relativistic field theory of two self-interacting Dirac fermions coupled 

to a Chern-Simons gauge field. The field theory has a phase diagram with 

the correct number of Goldstone modes in each regime and a phase transition 

at a critical coupling A* > Ac. We identify this transition with the isotropic 

Heisenberg point. It has a non-vanishing Neel order parameter, which drops 

to zero discontinuously for A < A*. 
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I. INTRODUCTION 

Since the discovery of high-Tc superconductors |l[ , the two-dimensional quantum Heisen- 
berg model has received considerable attention. This is largely due to well established ex- 
perimental facts which strongly suggest that these compounds can be described by a doped 
Heisenberg spin-1/2 quantum antiferromagnets 0. 

Dimensionality plays a crucial role in the properties of the Quantum Heisenberg Anti- 
ferromagnet. The 5 = 1/2 quantum antiferromagnetic chain can be solved exactly using 
the Bethe ansatz [^. By using the Wigner- Jordan transformation |^^, this model can be 
mapped onto a system of spinless, interacting, fermions with a coupling constant equal to 
the anisotropy parameter. This model is particularly simple in the XY limit where the spin 
problem maps to free fermions. Although fairly reliable in general dimensions, in one space 
dimension the spin wave theory is plagued by a number of notorious problems. This approx- 
imation is based on the Holstein-Primakov ^ transformation which maps S = 1/2 spins 
into hard core bosons. The spin-wave approximation ||^ relaxes the hard core constraint 
and treats correctly the commutation relation between spins in different sites. For one- 
dimensional systems, spin-wave theory (or rather, the 1/S expansion) is infrared divergent 
order-by-order. This divergence is the manifestation of the fact that the continuous sym- 
metry of global spin rotations cannot be broken in one space dimension. It also misses the 
essential fact that half-integer spin systems are critical while integer spin systems are always 
quantum disordered and have an energy gap [H . These properties of the exact ground state 
of the system can be recovered in one-dimensional spin systems by using non-perturbative 
methods, such as the Wigner- Jordan transformation combined with bosonization [p|,p!0|. 

Much less is known for two-dimensional quantum antiferromagnets. Firstly, there is 
no exact solution available in any limit of the spin-| system. Spin- wave theory predicts a 
Neel ordered ground state for the isotropic antiferromagnet on a square lattice, although 
with a moment reduced to 50% of the classical value by quantum fluctuations [|ri|]. The 
Hamiltonian for the anisotropic quantum Heisenberg antiferromagnet on a square lattice is 



H = J Y. |a5.(x)5,(x') + ^[5+(x)S-(x') + 5-(x')5+(xO]J , (1.1) 

<x,x'> '^ ^ '' 

where < x >,x' > denotes nearest neighboring sites and A is the anisotropy parameter 
(A = 1 corresponds to the isotropic case). 

The following facts are known to be true for this system. For A ^ 1, the Ising term 
dominates and the ground state should be close to a classical antiferromagnet which has 
total Sz = 0. This state has an energy gap and an expansion in powers of 1/A is rapidly 



convergent. In the opposite XY limit, where A — *> 0, there is a theorem ||T2[ due to Kennedy, 
Lieb and Shastry and to Kubo and Kishi, which proves that there exists long range order 
with the spins lined up on the XY plane for A < Ai (with Ai > 0.13). The same theorem 
proves that, in the Ising regime, the antiferromagnetic ground state extends at least down 
to an anisotropy parameter A < 1.78. No theorem is known for the isotropic case A = 1 and 
spin S = ^. For S > 1 Dyson, Lieb and Simon |T^ proved a theorem which shows that there 



is Neel order even at the isotropic antiferromagnetic point. Finite size diagonalization [IT^ 



quantum Monte Carlo |T^ and variational estimates |]TB[, are more consistent with a Neel 
antiferromagnetic ground state for the isotropic antiferromagnet. For two-dimensional quan- 
tum antiferromagnets, the semiclassical 1/5* expansion is free of the infrared divergencies 
found in one dimension. This approach predicts that the low energy limit of the isotropic 



antiferromagnet is a non-linear sigma model without a topological term |T^. This latter 



results have been confirmed by detailed renormalization group studies |jT^ which yield an 
excellent agreement with experiments on La2Cu04. 

These results suggest that the anisotropic quantum antiferromagnet has a phase diagram 
with just two phases: (a) a Neel state with Ising anisotropy for A > 1, and (b) an XY 
phase for A < 1. For A > 1 the Ising anisotropy should make all excitations massive {i.e. 
no Goldstone bosons for A > 1). For A < 1 the U{1) XY symmetry is spontaneously 
broken and there should be one Goldstone boson (spin wave). In this scenario, exactly at 
the isotropic point A = 1, the SU{2)/U{1) global symmetry of the Heisenberg model is 
spontaneously broken and there should be two Goldstone bosons (spin waves), as predicted 



by the non-linear sigma model [0. In some sense there is a phase transition at A = 1 in 
that the Neel order parameter should jump dis continuously to zero as A is decreased through 
A = 1. Precisely at X = 1 the XY and Neel orders are equivalent under an SU{2) rotation 
and, thus, there is still long range order. In contrast, the one-dimensional spin-| chain has 
a line of fixed points for A < 1 and Neel order in the massive phase A > 1. In section |I| we 
give an argument, based on the 1/5 expansion, in support of this general scenario. 

In this paper we investigate the anisotropic quantum Heisenberg antiferromagnet on a 
square lattice using a generalized Wigner- Jordan transformation constructed earlier by one 



of us [|T9|. In reference [|T9| it was shown that the Wigner- Jordan transformation in two 
dimensions is a special form of a statistics changing transformation which, quite generally, 
is achieved by coupling particles to (lattice) Chern-Simons gauge fields with a properly 
chosen coupling constant 6. The quantum Heisenberg antiferromagnet becomes equivalent 
to a system of spinless fermions which interact with each other (just as in the case of the 
spin chain) but which are also coupled to the Chern-Simons gauge field. Thus, even in the 
XY limit this system is interacting. Systems of this type have been considered recently in 



connection with the problem of anyon superfluidity [p0| -|2^ and the Fractional Quantum Hall 
Effect (FQHE) |2^J2^. Unlike the case of the spin chain, the equivalent fermion problem 
is never free and a new type of approximation has to be found. The analog of mean field 
theory in this context is known as the Average Field Approximation (AFA) (defined below 
in Section |T|). We will also show that the physics of this problem is hidden at the mean- field 
level and it is only revealed by a careful consideration of the fluctuations. 

Just as in systems of anyons or in the FQHE, the energy spectrum predicted by the AFA 
consists of free fermions with an effective band structure generated by the self consistent 
flux. We flnd that, for the sector with 3^ = (which corresponds to half-flUing of either 
the hard-core bosons or the fermions) the average uniform effective flux is equal to one- 
half of the flux quantum. Thus, on average, we flnd a flux phase analogous to that of 
the gauge theory approach to quantum antiferromagnets of Baskaran and Anderson ||25[| , 



Affleck and Marston ||2^ and Kotliar [^. The AFA also predicts that, for A > Ac (where 



Ac ~ 0.39), a gap opens up in the energy spectrum. In this regime the fermion density and 
the flux acquire a modulation with wave vector (7r,7r). The Wigner- Jordan transformation 
maps the 2;-component of the spin Sz{x) at x to the fermion occupation number n{x) by 
Sz{x) = I — n{x). Thus, we identify this regime with Neel antiferromagnetic order. For 



A < Ac the AFA spectrum of fermions is massless. Recently, Wang [Q has studied the 
Heisenberg model on a square lattice using the Wigner- Jordan transformation of reference 
||T9| combined with an approximation similar to the average field approximation discussed 



in Section |^ The results of the average field approximation (AFA) that we present here 
disagree with Wang's, mainly because his form of the AFA is not fully self-consistent. 

This mean-field spectrum is incompatible with the scenario proposed above, based on 
the semiclassical expansion. There, the spectrum of low lying states contains only integer 
spin fluctuations {i.e. spin flips). Some of these states may be massless ( as for A < 1) or 
all massive (such as for A > 1) but they are all bosonic. In particular, and in contrast with 
one-dimensional systems, the 1/S expansion predicts the existence of long range order. 

This problem is solved by a careful consideration of the role of the fluctuations. As 
expected, symmetry plays a crucial role here. The main problem is that, for small values of 
A, the AFA fermion spectrum is gapless. For arbitrary values of the Chern-Simons coupling 
6, two types of gaps, even or odd under Time Reversal (T) or Parity (P) can be generated 
by fluctuations. The Lagrangian at the level of the AFA is even under both P and T. Thus, 
we should expect that fluctuations will generate all terms with low scaling dimension [i.e. 
operators which are either relevant or marginal) which are compatible with the symmetries 
of the full system. Notice that here we encounter the opposite of the situation usually found 
with spontaneously broken symmetries where the mean-field theory has less symmetry than 
the full system. 

The symmetry analysis becomes more transparent in terms of an effective theory for 
the low energy degrees of freedom. This effective continuum theory for the two-dimensional 
quantum antiferromagnet on a square lattice turns out to be a theory of two species of 
relativistic Dirac fermions (moving at the Fermi velocity defined in Section 0) coupled to a 



Chern-Simons gauge field and to the fluctuations of an effective Neel order parameter field. 
This effective field theory is derived in Section |V[ This theory is a generalization to 2D of 
the well known equivalence between the antiferromagnetic spin chain and a field theory of 
interacting relativistic fermions in 1 + 1 dimensions, known as the Luttinger-Thirring model. 

In this language, the fermion spectrum can be understood in terms of the possible energy 
gaps, or masses, of the fermions and of the symmetry properties of the operators connected 
with these masses. Our analysis shows that the Neel order parameter acts like a mass 
operator which does not break P or T and which we will denote as MnsgI- In the Neel phase 
the two species of fermions acquire a mass, but with relative opposite sign. However, for 
general values of the coupling constant 9, the Chern-Simons term breaks both T and P. 
Thus, quantum fluctuations of the Chern-Simons field will necessarily generate all terms 
which break the same symmetries. A fermion mass term with the same sign for both species 
also breaks both P and T and we find that it does get generated by quantum fluctuations. 
We will refer to this as the induced fermion mass. Mind- The actual phase diagram follows 
from the competition of these two mechanisms. 

The generation of a parity breaking mass term changes radically the long distance be- 
havior of the system. It is well known pS that fermions with masses which break time 



reversal and parity induce Chern-Simons terms in the action of the gauge fields at length 
scales long compared with the correlation length of the fermions, i.e. the inverse of the 
parity breaking fermion mass. Also, the sign of the induced Chern-Simons term is equal to 
the sign of the mass of the fermion. Hence, at length scales long compared with 1/Mind, we 
expect to see a finite renormalization of the Chern-Simons coupling constant from its bare 
value 9 to some effective value. This effective value depends on the pattern of symmetry 
breaking, i.e. on the relative signs of the induced masses which, in turn, are determined by 
the bare coupling 9 itself and by the nature of the ground state. We find that the sign of 
Mind is such that the induced Chern-Simons coupling tends to reduce the bare Chern-Simons 
coupling. This renormalization can be viewed as a tendency to screen the bare statistics. 
This is a manifestation of a more general property which we may think of as a "Lenz Law 
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of Statistics" iSOl. 



Thus, we arrive to the following scenario. At small A, the fermions acquire a parity 
breaking mass Mjnd through the fluctuations of the gauge field. In turn, at distances long 
compared with l/Mjnd, a Chern-Simons term is induced with a coupling constant which tends 
to cancel the bare statistics. We will find in Section [V^ that for 9 = Xjlii the cancellation 
is complete and the gauge fields are actually gapless at long distances! This scenario is 
reminiscent of anyon superfluidity. We will identify the gapless transverse gauge fluctuation 
[i.e. Laughlin's mode ^0|) with the gapless transverse spin wave of the XY regime. Clearly 



two fermion masses with different symmetry necessarily compete with each other. We expect 
that when they become of comparable magnitude |Mind| ^ |MNeei| a phase transition should 
occur. For this range of A, one of the two species of fermions has a mass that is very small 
and, at least nominally, is going to vanish at some critical value of the anisotropy A*. We 
will identify this phase transition with the S'f/(2)-symmetric Heisenberg point |^T|. This 
phase transition occurs at a value of the anisotropy A* > Ac. This phase transition preempts 
the naive second order transition predicted by the AFA from taking place. This mechanism 
seems to bear a close analogy with a fluctuation-induced first order transition P^J53[ . 



To summarize, the Chern-Simons (or Wigner- Jordan) approach to the spin S = ^ 
anisotropic quantum Ifeisenberg antiferromagnet yields a phase diagram which is in qual- 
itative agreement with the predictions of the semiclassical 1/5* expansion. The physically 
correct (P-C.) phase diagram is not found at the level of the Average Field Approximation 
and, except in the regime of strong Ising anisotropy, it is due almost entirely to fluctuation 
effects. In this paper we show how this physical picture is realized in the context of the 
Chern-Simons theory. This is a (necessary) first step before these methods could be applied 
to more subtle problem such as frustrated systems. 

This paper is organized as follows. In Section |1| we develop the semiclassical 1/S theory 
of the anisotropic antiferromagnet. In Section JTT| we review the Wigner- Jordan construc- 
tion and its connection with the Chern-Simons theory on a square lattice. In Section |rV| we 
present the results in the Average Field Approximation and discuss the role of gaussian flue- 



tuations. In Section the effective field theory is derived. The dynamical, non-perturbative, 
effects of ffuctuations are discussed in Section [V^ where we give a justification of the phase 
diagram discussed in this (long!) introduction. Section |Vll| is devoted to the conclusions. 



II. THE ANISOTROPIC HEISENBERG QUANTUM ANTIFERROMAGNET 

FOR LARGE S 

In this Section we discuss the semiclassical 1/5* theory of the anisotropic antiferromagnet 
in a square lattice. 

The easiest way to get a path integral quantization for a spin system is to use coherent 
states. In this Section we will follow the methods described in reference [|5 . 



The set of coherent states {\n >}, labelled by the unit vector n, is generated by a rotation 
of the highest weight vector (IS*, S >) of an irreducible representation of the group SU(2) of 
spin s of the form 

|^>=gie(noXn).5|^^^^ (2.1) 

where Uq is a unit vector along the quantization axis, 6 is the colatitude ( n.no = cos^) 
and Si {i = 1, 2, 3) are the three generators of SU(2) in the spin-s representation. The state 
\n) can be expanded in a complete basis of the spin-s irreducible representation {\S,M >} 
where M labels the eigenvalues of S3. The coefficients of the expansion are the representation 

matrices D^^\n)M,s 

s 
\n>= Y. D^'^\n)M,s\S,M> (2.2) 

M=-S 

The matrices D^^' do not form a group but satisfy 

D(^)(rii) D(^)(n2) = /^(^^(ns) e'*("i."2,-3)53 (2.3) 

where ni, n2, and ns are three arbitrary unit vectors and $(f2i, 722, ^^3) is the area of the 
spherical triangle with vertices at ni, n2, and n^. Other useful properties of the spin coherent 
states are: (a) the inner product < ni\n2 > 
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(b) the diagonal matrix elements of the generators S 

< n\S\n >= sn (2.5) 

and (c) the resolution of the identity operator 

d^{n) \n >< n\ (2.6) 



where we have used the integration measure 



2s + l\ ,3^. ,^2 



dij{n) = ( ^^ J d-'n6{n' - 1). (2.7) 

Using these properties we can write an expression for the path-integral in this coherent state 
representation. The zero temperature partition function reads 

Z= /pne^'^"!"] , (2.8) 

where the action for the many-spin system in real-time is given by 

/•T 

<SM[n] = s^Swz[n{r)] - I dx^J-s^ ^ {n^{f,XQ).n^{r' ,xq) + An3(f,xo).n3(r",xo)} 

r <r,f'> 

(2.9) 

To write this expression we have used the Hamiltonian for the anisotropic quantum Heisen- 
berg antiferromagnet on a square lattice given by eq. (|1 . 1|) (with < r,r' > denoting nearest 
neighboring sites), and we have assumed periodic boundary conditions. In Eq. (|2.9|) n±_ is 
the projection of the vector n onto the 12 (or xy) plane. 

The first term in eq. ( p.9|) is just the sum of the Wess-Zumino terms (or Berry phases) 
for the individual spins. The contribution of each term to the action is 

Swz[n]= / dr / dt n{t,T).{dtn{t,T)xdrn{t,T)) (2.10) 

JO Jo 

where jS = iT. The effect of this term is to quantize the spin. 



The effective action iSa/ [ti] scales like s, the spin representation. Therefore, in the large 
spin limit, the path integral eq. ( p.8[ ) should be dominated by the stationary points of the 
action. This is the semiclassical limit. Corrections to the large-s limit can be arranged in 
an expansion in powers of -. Since we expect to be close to a Neel state, we will stagger the 
configuration 



n(r] 



-l)^i+^2n(r) (2.11) 



The Wess-Zumino terms are odd under the replacement of eq. (|2.11|) and thus get staggered. 
Up to an additive constant the action reads 

f 

+A|[n3(f,xo) -n^{f+ej,XQ)f + [n^{f,XQ) - n3(f- e^, xq)]^} 

+ (l-A){2[n3(f,xo)]' + [n3(f+e„xo)]' + [r?3(r-e„a;o)]'} (2.12) 

We split the staggered spin field n in the following way 

n{r) = m{f) + {-iy'+''^aol'{f) (2.13) 

where rh{r) is a slowly varying piece , the order parameter field, and l{r) is a small rapidly 
varying part which roughly represents the average spin. The constraint fp = 1 and the 
requirement that the order parameter field fh should obey the same constraint, m^ = 1, 
demand that rh.l = 0. Using this property we can write the lagrangian density for this 
theory in the long wavelength limit as 

CMi'niJ) = — l.{mxdo'm) 
«o 

-Js' I Y: mrn^)' + Hd,m,r] + 8[{Lr + XihY] + 4(1 - A)[^ + {^1 (2.14) 

[j=l,2 «0 J 

where oq is the lattice spacing. 

In the long wavelength limit, the Wess-Zumino action can be written as a sum of a 
topological term and the first term in eq. ( p.l4|) . It has been shown ( |T^) that if we expect 
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to have Neel order, the topological term does not contribute to the action in two space 
dimensions. Notice that, in the one-dimensional case, this same procedure leads to a sigma 
model with a topological term. 

After integrating out the fast modes, i.e., the components l± and I3 of I, in the partition 
function, the resulting lagrangian is 

CM{rn) = — [— (<9om±)^ - Vs Y] {djm±f] 

^^^ " ^^ ^^ im^r + ^^-^^{m.ndom^r - {dom.r] (2.15) 



al 2g'-~^' ■ 2 2gv 
where the coupling constant g and the spin wave velocity Vg are given by 

g = -ao{l + \)-^ (2.16) 

s 

v, = 4Jsao(l + A)5 (2.17) 

We can Wick-rotate back to imaginary time (i.e. X3 = ixo), and write the Euclidean La- 
grangian density Ce as 

^g Vg j=l,2 "^9 ^'^s i=l,2 

■ ^^^ " ^^ ^^^^3)^ + ^^^^{rTisfmrn^r - idsms)'] (2.18) 



al 2g" '"' ' 2 2gvs 
By direct inspection of eq. (|2.18| ) one can see that the third term of this action is relevant in 



the long wavelength limit. The physics of this term is the following. For A < 1, the system 
will lower its energy by making rris -^ 0, i.e., the system is in the XY limit. On the other 
hand, if A > 1 the energy will be maximized when m-^ acquires its maximum possible value, 
i.e., m^ = 1 and the system is in the Neel state with Ising anisotropy. 

The phase diagram for the anisotropic quantum antiferromagnet suggested by these 
results has only two phases. For A > 1 the system is in a Neel state with Ising anisotropy, 
and all the excitations are massive. For A < 1 the systems is in an XY phase, the U(l) XY 
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symmetry is spontaneously broken and there should be one Goldstone boson. Exactly at 
the isotropic point (A = 1) there should be two Goldstone bosons as predicted by the non 
linear sigma model |T7| . 

III. CHERN-SIMONS ON A LATTICE 

We begin by reviewing the path integral picture of a spin system on a two dimensional 
lattice in terms of fermions coupled to Chern-Simons gauge fields, introduced in reference 

El- 

The Wigner- Jordan transformation is based on the identification of a system of hard core 
bosons {i.e. spin flips) with an equivalent system oi fermions each of them rigidly attached 
with solenoids that carry one-half of the flux quantum. Mathematically, the equivalent 
system is a theory of fermions coupled to a Chern-Simons gauge field on the square lattice. 



Chern-Simons theories [35| have been used with great success in the Fractional Quantum 



Hall Effect p3| , p^ and in anyon superfluidity [p0| -p2| . The presence of the lattice introduces 
a number of subtleties not found in continuum systems. The role of the Chern-Simons 
gauge fields is to enforce the constraint that attaches particles to fiuxes locally, and a set of 
commutation relations among the gauge fields compatible with these constraints |^. These 
two features are key ingredients for the Wigner- Jordan transformation to work. 

However, unlike the one-dimensional spin chain, in two dimensions the equivalent system 
of fermions is coupled to a gauge field which can have local fiux. Hence the fermions are 
always interacting, even in the XY limit, and approximations become necessary. Written 
in the fermion language, the system can then be described in terms of a theory of fermions 
which interact with each other and with a Chern-Simons gauge field. 

In what follows we will use a path-integral description. The zero temperature partition 
function for this problem has the form 

Z= fvilj*ViljVAf,e'^ , (3.1) 
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where %l){x,t) is a fermi field {i.e. Grassmann variables in the path integral) defined on the 
sites {x} of the square lattice and A^ are the statistical or Chern-Simons gauge fields. The 
space components Aj{x,t) are defined on the links of the lattice while the time component 
Ao{x,t) is defined on the sites. The role of the gauge field is to change the statistics. The 
action S is given by 

+2A(|^(x,t)r-i)(|V^(x + e„t)r-i))} + ^^,3M) (3.2) 

where Dq = d{) + iAo is the covariant time derivative, and the spatial covariant derivative is 
in this case the gauge covariant lattice difference implied by the hopping term. The action 
S of Eq. ( |3.2| ) describes self-interacting fermions which are coupled to a fluctuating Chern- 
Simons gauge fleld. The self-interaction is represented by the quartic term in fermions in 
the action and it corresponds to the SzSz Ising interaction of the Heisenberg model. We will 
call this term Si^t. 

The lattice Chern-Simons Scs {A) action was deflned in references p!9|j36|] . Its explicit 



form will be given below. The coupling constant 6 is chosen to be 6' = 2~: so that the 
statistics corresponds to bosons (with hard cores). For general Chern-Simons coupling 9 
this system is equivalent to a system of interacting anyons with statistical angle 5 = ^, on 
a square lattice |T^. Lattice anyons have been studied numerically by Canright et. al. ||3^ 



and analytically by Fradkin |2^ . We will see in Section |^ that the problem at hand is an 



example of the degenerate solution found in reference p2 . 

In the representation of the Heisenberg model in terms of fermions coupled to gauge 
fields, with the action of eq. ( |3.2| ), the natural mean field approximation consists of detach- 
ing the fermions from their local fiuxes and to replace this dynamical fiux by a static average 
background. Unlike spin-wave theory, in this mean field approximation the hard core con- 
straint is taken into account exactly. The phase factors present in the hopping amplitudes of 
the equivalent system of fermions, whose role is to enforce the original bosonic commutation 
relations, are treated approximately. The Average Field Approximation (AFA), as this mean 
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field theory has come to be known, was first introduced by Laughlin in the context of his 
study of the anyon gas [^ and subsequently used quite extensively in the context of the 
FQHE by two of us [^ . A peculiar feature of this mean field theory is that it breaks a num- 
ber of space-time symmetries in a very explicit manner. For example, in the anyon gas the 
ground state obtained at the AFA level breaks Galilean invariance while the actual ground 
state does not. The fiuctuations around the AFA restore Galilean invariance. Likewise, in 
the context of the FQHE Galilean (or rather, magnetic) invariance is broken at the level of 
the AFA but it is also restored by gaussian fiuctuations [^. 

As stated above, the role of the Chern-Simons gauge fields is to change the statistics from 
(hard-core) bosons to fermions. Here we will follow the approach of references [0,^]. The 



effect of the Chern-Simons action is twofold: (a) a constraint on the allowed states which are 
required to satisfy a relation between the local particle density and the local statistical fiux 
and (b) a set of commutation relations for the gauge fields |T^. With the sign conventions 



of the action of eq. (p.2|), the constraint reads 

p{x,t) = eB{x,t) (3.3) 

with p{x, t) = ip*{x, t)ip{x, t). This is a constraint on the allowed states in the Hilbert space 
and it plays the same role as Gauss' Law in Maxwell's electrodynamics. In eq. ( p.3|) the 
particles live on the sites of the square lattice whereas the fiux B{x, t) is defined on the sites 
of the dual of the square lattice , i. e. the center of the plaquette "north-east" of the lattice 
site X. The Chern-Simons gauge fields are defined to be on the links of the square lattice. 

The lattice form of the Chern-Simons action Scs (A) can be written as the sum of two 
terms 

-^cs (A) = S^ + 5(f (3.4) 

where S^l' and S^^' are responsible for enforcing the constraint and for the determination 
of the commutation relations respectively. Here we will use the form given by Eliezer and 
Semenoff ^6\ 
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"^cs^ = j dtY^Ao (f , t) e'^diAj {x, t) (3.5) 

X 

^i? = -\j dtY^A, (f , t) K^^^^A, (f , t) (3.6) 

X 

Here, we have used the forward difference operator rfj which acts on functions f{x) defined 
on the sites as dif{x) = f{x + ii) — f{x), where Cj is a unit vector pointing towards the 
direction z = 1, 2 of the square lattice . Similarly the backward difference operator di acts 
like di = f{x) — f{x — ii). The kernel K^^ is found to be given by the matrix [^ 



K'^ = ~l 



I ^„ I ^_ O O^. I o^_ I ^_^. \ 



^2 + c?2 -2 - 2di + 2d2 + d2d 



\ 2 + 2d2 — 2di — did2 ~di — di i 



(3.7) 



The quartic term in the fermion part of the action represents the SzSz interaction. The 
constraint of Eq. ( |3.3|) restricts the space of configurations to those in which the fermion 
occupation number at a site is equal to the flux at the plaquette north-east of the site divided 
by 9. Hence, it is legitimate to replace in the action the fermion density by B /9. Therefore, 
the interaction term in the action, S'int, becomes only a function of the configuration of the 
gauge fields 

5int = -^E(^^(x,t) - ^)V(x-x')(^0(x',t) - \) (3.8) 

x,x' 

with a pair potential ^(x — x') given by 

\ J\ if x' = X ± e, (j = 1, 2) 

r(x - x') = . u , ; ^3 g^ 

otherwise 

Notice that the interaction term is bilinear in the gauge fields instead of a quartic functional 
of the Fermi fields. This result was obtained before, in the context of the FQHE, in reference 
Alternatively, one could use a Hubbard-Stratonovich transformation and arrive to the 



same result 38 



By putting all the terms together we arrive to the final form of the action 

s = s^iip.r.A,) + SUA) + ds^iA) (3-10) 
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where S'f('?/','?/'*,^^) is the action for the fermions coupled to the gauge field 

-'-^^■("•*)'*V(x,t) + c.c. 






(3.11) 



Smti-^fi) and Scs are defined in Eqs. ( |3.8| ) and ( |3.4| ) respectively. From now on we will use 
the action of eq. ( p.lO|) . 



IV. MEAN FIELD THEORY AND SEMICLASSICAL EXPANSION 

We will now proceed to derive a mean field theory in the usual fashion. The fermionic part 
in the action (|3.1CI|) , being bilinear, can be integrated out yielding a fermion determinant. 
The resulting effective action ^efr is given by 

5eff = -itT\og[iDo + /i - h{A)] + eScsiAf,) 

-l/rftE(^S(x,t)-i)V(x-x')(0S(x',t)-i) (4.1) 

x,x' 

where h{A) is the kinetic part that we can write in operator form as 

M^) = ^EE |x,t)e*'^^('^'*)(x + e„t|+h.c. (4.2) 

^ X j=l,2 

The semiclassical approximation is obtained by expanding around stationary configura- 
tions of fields that minimize the effective action. It is worthwhile to note that this effective 
action does not contain any small parameter to control this expansion. This is a problem 
that was also found in the context of the anyon superfiuid as well as in the fermion Chern- 
Simons approach to the FQHE of references p^j38|] . There it was found that if the AFA 
had a spectrum which is fully gapped, the fiuctuations restore the correct spectrum at long 
wavelengths. We will see that, for the problem at hand, the AFA yields a gapless spectrum 
at least for a range of values of the anisotropy parameter A. Thus, the validity of the AFA, 
even qualitatively, is questionable for that range. Indeed we will find it necessary to go well 
beyond the AFA in order to obtain asymptotically exact results. 
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The Average Field Approximation is realized by the solutions of the saddle-point equa- 
tions 

t4¥^ - = (4.3) 

As usual, the variations of the fermionic part S-p of the action S^s {i.e. the first term in 



eq. ( [4. ID ) with respect to the components of the gauge field Afj, gives AFA expressions for 



the charge density ?2(x, t) 



(n(x, t)) = (^*(x, t)^(x, t)) = --o^ (4.4) 



and current density jki^,t) 



,%J r 



Ofc(x,t)) = (^ [^*(x,t)e'-^'=(^'*V(x + efc,t) -^*(x + efc,t)e-^-^'=('''*)^(x,t)J) 



Within the AFA, we find that the average density and currents are given by 



5S- 



F 



5A(x,t) 

(4.5) 



(n(x, t))AFA = - . . / ,^ = -«5(x, t; X, t) (4.6) 

Mo(x,t) 

(J.(x, t))AFA = -J^^^ = ^ {S{^ + e,, t; X, t)e^^^(-'*) - ^(x, t; xe,, t)e-^-^^ (-'*)) (4.7) 
where Dq = do + iAo and ^^ (yU = 0, 1, 2) is the expectation values of the components of the 



gauge fields within the AFA. The function S{x,t;x.',t'), which appears in ( ^I6D and ( [4.7|) , is 
the Green function for the fermions moving in the background field Afj, which is the solution 
of the lattice differential equation 

tDo + l2- h{A)) ^(x, t; x, t') = (5x,x'<5(t - t') (4.8) 

Below we give the explicit form of this Green function. 

By varying ^efr with respect to ^o "we recover the constraint equation, now as a condition 
for the value of the local density of the stationary configurations 

(n(x)) = e{Bi^)) (4.9) 
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Likewise, by differentiating with respect to the spacial components Ak, we find an equation 
for the fermion current 

Ofe(x)) = 9e''d,{Ao{^)) + I [k'' - K"') 9o(A(x)) - Oh^i E V{^ - x')rfKS(x')) (4.10) 



In terms of the pair potential of Eq. ( |3.9| ), we can write the expectation value of the current 
in the form 



Ofc(x)) = 9e'^di{Aoi^)) + - (k'' - K'') d,{Ai{^)) 



e'JXekidi Yl ((^(x + e,)) + (i3(x - e,))) 
i=i,2 



(4.11) 



In Eqs.. ( 4.10 ) and ( 4.11 ) K^'' is the operator matrix defined in Eq. (|3.7| ). Using this definition 
explicitly, the terms in Eq. ( 4.11|) which depend on i^'^ are given by 



e[K"'-K''')do{Ai{^)) = --€''' 



4 + 2(^2 - 4) + 2(di - di) + d2di + rfidzj do{Ai{^)) 

(4.12) 



In the continuum limit, the first two terms in Eq. ( [4.11|) are equal to the conventional 
Chern-Simons current, i.e. 6e'''-{Si{'x)). 

The Saddle Point Equations ( [4.9|) and ( [4.11|) have many solutions. For a half-filled 
system {i.e. 3^ = 0), the solution with lowest energy corresponds to a stationary state with 
a modulated charge density and with zero current. In the magnetic language this is a state 
with a non zero Neel order parameter A, such that 



(n(x)) = - - Ae^-\ 



(4.13) 



Using the constraint of Eq. ( ^4.9| ) we get 



(S(x)) = -(n(x)) 
= 27r(n(x)) 



(4.14) 



This state is time independent and it does not support any current, local or global. For a 
square lattice with N x N sites (with N even), periodic boundary conditions and 6 = 1/2tt, 
we can satisfy this static constraint with the following choice of gauge fields 



Ao = ioe'"-^ A(x) = Sj,2 (f + vf • f + TrAe^'^-^) (4.15) 

This solution corresponds to a problem in which a fermion moves in the presence of a 
modulated magnetic field (with an average of half flux quantum per plaquette) and a periodic 
potential A, with the same modulation. By solving these equations we find a solution in 
which 

Ao = 4AJA (4.16) 

For ground states with a modulation in the effective magnetic field, a periodic site potential 
(^o) appears which is commensurate with the variation of the field. The charge and the 
field vary in the same way, as required by the constraint. Thus, everything is self-consistent. 
Let us note in passing that for sectors with Sz y^ the average fiux is not equal to one-half 
of the quantum. At the level of the AFA this problem now becomes equivalent to a general 
Hofstadter problem [39|. This problem is known to a have a very complex spectrum which 



we will not explore here. One of us |^2| discussed a similar problem in his treatment of the 
lattice anyon gas. In what follows we will only discuss the case 5*^ = 0. 

The Green's functions of the saddle point problem are obtained by solving the lattice 
differential equation of Eq. (^4.8| ). In the gauge that we have chosen, ^i = and with the 
configuration of gauge fields of Eq. (^4.15|) , the solution of Eq. ( [4.8|) is the Green function of 



a one-particle system with the effective Hamiltonian Hmf = h{Aj) + X)x-^o(x), which takes 
the form 

^MF = E I^A^^e'"" |x)(x| + ^ (|x)(x + ei|+a(x)|x)(x + e2|+h.c.)| (4.17) 

with 

a(x) = e^-^^W = -sinvrA + ie^'^-''cos7rA . (4.18) 

Note that a(x) is equal to a = ie^'^^ when x belongs to the A sublattice and to a* = — ze~"^^ 
when X belongs to the B sublattice. For A = we recover the flux-phase state of (uniform) 
half flux quantum per plaquette pGJ^ . 



19 



The Hainiltonian of Eq. ( ^4.17|) can be diagonalized in Fourier space. The charge density 
and the effective flux are periodic functions which take different values on the sublattices 
A and B. Thus, the Fermion Green function is a matrix whose entries label the sublattice 
dependence of its arguments. It has the form 



f^apyXj X ) 



Saa{x,x') Sab{x,x' 
Sba{x,x') Sbb{x,x' 

where x = (x, t) and a, /3 = A, B. We find 

/ 



(4.19) 



Sap{x,x') = I 



^°° du 



dk e«^(*-<')~«(x-x') 



oo 271 J \ki\<^ tt'^ uj'^ - E'^Ck) +ie 



u + 4 JAA J(cos ki + a cos k2 

J {cos ki + a* cos k2) uj — AJXA 

(4.20) 



Notice that the momentum integrals are done over the range |fcj| < | (z = 1, 2). The fermion 
dispersion E(k) is given by 



E{k) = + JJcos^ ki + cos^ k2-2 sin(7rA) cos ki cos ^2 + (4AA)^ 



(4.21) 



Recently, Wang [^ obtained similar results but without self-consistency between charge 
and flux modulations. 

The dependence of E{k) in ttA as given by eq. ( 4.21|) is a consequence of the self- 
consistency. Once the Green function is known, the parameter A can be calculated by 
demanding that the saddle-point equations be satisfied. Thus, the average density on one 
sublattice, say A, must be given by 



/+00 f^^ r 
-oo Zn Jn 



dk UJ + AJXA 

|fci|<f TT^ uj"^ — E'^{k) + ie 



1 ^,, ^ f dk 1 

--2JAA / ^T^TTT 

2 i|fc,|<f 7r2 E{k) 



From this expression we arrive to the gap equation 

, . . dk 1 

A = 2JXA 



(4.22) 



fc,|<f 7r2 E{k) 



(4.23) 
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We now discuss the properties of the spectrum found in the AFA. For general values of A, 



the spectrum of Eq. (|4.21|) has a gap Eg = 4JAA at the four points in k space (±7r/2, ±7r/2). 



Eq. ( [4.13| ) is a self consistent equation for A. Qualitatively, its solution as a function of the 
anisotropy parameter A is a monotonically increasing function which begins at a critical 
value of A. From Eq. ( ^.231 ), it follows that there exists a critical value of the anisotropy 



parameter Ac given by 



1 r dk 1 



_: / — , ~ 2 X 1.285 = 2.570 , 4.24 

Ac Ak.\<^ T^' Jcoskj + cosk^ 



below which A = 0, and the spectrum is gapless. Thus, the AFA predicts that the spectrum 
of fermions has a gap above a critical anisotropy Ac — 0.39. For small and positive values of 
A — Ac, the gap has the dependence 

A~const.(A-Ac) (4.25) 

and it vanishes for all A < Ac. 

The vanishing of the gap for A — ^ Ac, and in particular the exponent of eq. ( [4. 25] ), results 



from the collapse of the Fermi surface to four Fermi points and from the linear dispersion 
near the Fermi points. Since the density of one-particle states vanishes in the middle of the 
band {E = 0) for a system with a linear, relativistic-like, energy-momentum dispersion in 
two space dimensions, all instabilities are pushed to finite values of the coupling constants 
and there is a critical coupling. In contrast, in conventional mean field theories of interacting 
fermions at finite density on a lattice (but with zero gauge flux) there is no critical coupling 
and the spectrum is always gapped in the presence of nesting. In such cases, the gap 
dependence for small A would be of the form A ~ e~^ . Thus, the existence of the critical 
value Ac is a consequence of the inclusion of the gauge flux that removes the van Hove 
singularity in the density of states characteristic of the two dimensional square lattice. 

The exponent of eq. ( |4.25| ) is valid only at the level of the AFA. In critical phenomena, it 
is usually the case that the exponents found at the level of "classical" approximations such 
as mean field theory or in the large- iV limit, get modified due to the effects of fluctuations. 
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The AFA is a semiclassical approximation. Since the dimensionahty of space-time of this 



system is 2 + 1 we should expect non-classical behavior. In fact, the exponent of eq. ( f4.25| ) 
would be correct for a theory of A^ species of self-interacting relativistic fermions in the 
N ^ oo limit. We will see in section |V| that the system that we are studying here is indeed 
related to a theory of self interacting relativistic fermions but not in the large A^ limit. 
Thus, in principle, fluctuations are expected to correct this exponent. However, we will also 
show that this second order phase transition is never reached and that it is preempted by 
a fluctuation induced first order transition at a value of lambda A* strictly larger than Ac. 
We will also show that, at this fluctuation induced first order transition,the system has the 
expected physical properties of the isotropic Heisenberg antiferromagnet in a Neel state. 

When applied to the isotropic case (A = 1), the results of the AFA imply a Neel state 
with a value for the order parameter A = 0.442 and the energy E = — 0.314J per bond. 
These results should be compared with the best numerical estimate of the energy, — 0.334J 



per bond [0,^]. Nevertheless, it is interesting that this approximation yields a Neel state 
instead of a flux phase, as one might have guessed beforehand. 

But, is this the correct spectrum? According to the AFA, the spectrum consists of free 
fermions with an energy gap for A > Ac but gapless otherwise. Clearly this spectrum has 
nothing to do with what we found semiclassically in section |T|. Earlier work in anyons and in 
the FQHE suggests that fluctuations should play a crucial role. However, notice that since 
for A < Ac the AFA spectrum is gapless, the fluctuations may yield much more important 
effects than what we have described with the AFA. In particular, this gapless case was found 
in reference ||2^ in the case of semions at half filling and was found that the state was not 



determined by the AFA alone. 

Hence, the next logical step is to look at the effects of fluctuations around the solutions 
of the AFA. While it is possible to carry out this calculation and to keep the full lattice 
effects at the same time, the expressions are very cumbersome and not amenable to an 
analytic treatment. Instead, we will resort to a different approach in which only the low 
energy degrees of freedom are kept. This is equivalent to replace the system by an effective 
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continuum field theory. These methods are accurate provided that the gaps do not become 
too large, in which case the lattice effects may become dominant. This limitation will force 
us to work at values of A close to Ac. Nevertheless, our results will be qualitatively correct 
even away from this regime. This approach is pursued in section|V|. 

V. EFFECTIVE FIELD THEORY NEAR Ac 

In this section we will consider the behavior of the system near the critical anisotropy Ac- 
The mean field theory of section |^ yielded an energy gap for the fermions which vanishes 
linearly as A — * Ac from above. The problem that we want to address here is the nature 
of this phase transition. We will find that the transition at Ac does not take place as a 
result of radiative corrections, namely of fluctuations of the gauge field. Instead, a different 
transition with a larger critical value A* will actually occur and supersede the AFA-predicted 
transition at Ac. By counting Goldstone modes we will be able to identify the transition at 
A* with the isotropic Heisenberg point. 

We will now develop an effective theory for the low energy degrees of freedom and 
construct an effective field theory for the fluctuations about the AFA. To do so we will adopt 
a point of view which follows quite closely the treatment of the one-dimensional anisotropic 
spin-1/2 Heisenberg antiferromagnetic chain 0. However, the physics that we find is very 
different. The Jordan- Wigner approach that we have used in this paper closely resembles 
the Jordan- Wigner transformation of one-dimensional systems. The main difference is that, 
in one space dimension, the only flux that can be defined is the one that is trapped by 
the entire chain and, hence, it is equivalent to a boundary condition. In two dimensions, in 
addition to global or topological flux, it is possible to generate local flux. Local fluxes cannot 
be reduced just to boundary conditions and a local, fluctuating, gauge degree of freedom 
appears necessarily in the effective theory. 

We begin by first summarizing the standard procedure used in one-dimensional systems, 
first introduced by Luther |^ . We will next follow that construction for the two dimensional 
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case and use it to discuss the critical behavior near Ac- 

The fermionization of the one-dimensional spin chains is usually done in the following 
manner p. First, the algebra of spin one-half operators is realized in terms of canonical 
fermion operators via the use of the Jordan- Wigner transformation. The resulting system 
consists of a set of spinless fermions hopping between the nearest neighboring sites of a one 
dimensional chain of atoms. Two fermions cannot share the same site (Pauli principle) and 
only interact when on nearest neighboring sites with a coupling constant equal to twice the 
strength of the SzSz coupling constant. The total fermion number is equal to the number of 
down spins (depending on the definitions). There are some subtleties concerning boundary 
conditions which are important but are not related to the issues that are being discussed 
here. Thus, the sector of the spin system with total 3^ = maps onto the half-filled 
portion of the Hilbert space of the Fermi system. Notice that in one dimension the fermion 
degrees of freedom exhaust the Hilbert space and no other degrees of freedom are needed 
to represent the states of the spin chain. The second step consists of finding a quantum 
field theory in one space dimension which yields the exact long distance critical behavior of 
this system of interacting spinless fermions |1T[]. This is done by separating the slow from 



the fast components of the Fermi fields. In one dimension, the fermions can either move 
to the left (left movers) or to the right (right movers). The non- interacting Fermi system 
(equivalent to the spin one half chain) has two Fermi points with momenta p = ±pi? = ±|. 
The left and right moving components of the Fermi field can be thought of as the chiral 
components of a two-component Dirac spinor in one space and one time moving at a "speed 
of light" equal to the Fermi velocity vp- The equivalent quantum field theory has the left 
and right movers interacting through a backscattering process (up to umklapp processes 
which are crucial to reproduce the correct behavior of the isotropic Heisenberg model [p^). 
The resulting field theory is the well known Luttinger-Tomonaga-Thirring model |^ which, 
by using bosonization, can be shown to be equivalent to the sine-Gordon field theory 
(if the umklapp terms are kept). The result is that for weak Ising coupling {i.e. strong 
XY anisotropy) the spectrum of the system is that of fermions with anomalous dimensions 
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(as a result of the backscattering interactions) up to a value of the backscattering coupling 
constant at which the umklapp processes become marginal operators (and beyond which, 
in the Ising phase, they are relevant). In the Ising phase there is an energy gap for all 
excitations. 

In the case of the two-dimensional spin system, the AFA of section |V] yielded a spectrum 
of free fermions with a "flux phase" band structure. At half filling ,i.e. total Sz = 0, the 
flux phase band structure has a "Fermi surface" which reduces to four Fermi points located 
at (|, |) , (— f , f ) , (f , — f ) and (— f , — f )• Just as in the case of the one- dimensional spin 
chain, whose fermion description has two Fermi points, the two-dimensional problem can 



also be mapped onto a Dirac-like problem of Affleck and Marston [^]. The main differences 
between the problem that we discuss here and the flux phase of Affleck and Marston is that 
(a) the fermions here are spinless and (b) instead of a Hubbard-Stratonovich field which 
fiuctuates both in amplitude and phase we have just a fiuctuating phase on the bonds. The 
Chern-Simons term is not present in the system discussed by Affleck and Marston. 

We now follow the methods and notations of reference |^ to obtain an effective con- 
tinuum theory. Since the detailed derivation is rather tedious we will only highlight the 
procedure. Our general strategy will be to look for the terms in the action which involve 
only low energy degrees of freedom. We will only keep terms with the smallest numbers 
of derivatives in each of the ffuctuations since higher derivative terms are irrelevant in the 
renormalization group sense. We should keep in mind that this procedure will only give an 
approximate value for the coupling constants of the effective low energy theory since instead 
of integrating out the high energy degrees of freedom we are simply neglecting them. Hence, 
even though the form of the effective action will still be correct the values of the parameters 
[e.g. the anisotropy A) at which transitions may occur will not coincide with the predictions 
of the lattice system. In particular the critical anisotropy will not be precisely Ac = 1, 
but only approximately. We will have to use the symmetry properties of the spectrum of 
excitations at a certain value of the coupling constant of the effective theory to identify the 
isotropic point. 
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The starting point is the action of eq. ( P-lOj ). The flux phase that we found in Section |^ 
has a spectrum of fermions which become gapless at four Fermi points. The physically 
important states are those close to the Fermi points. Out of these states, we construct 
two (two-component) Dirac spinors. We begin by defining first a set of spinor components 
on the sites of the real square lattice. It is convenient to split the square lattice into four 
sublattices 1, 2, 3 and 4. Sublattice 1 is the set of sites of the form {x = (2ni, 2n2)} (with 
rii and ^2 arbitrary integers). Sublattices 2, 3 and 4 are the sites of the form {x + ei}, 
{x + 62} and {x + ei + 62} respectively. The fermion amplitudes on each sublattice will 
be denoted by ^a(x), with a = 1,...,4. Likewise, the gauge fields have to be split into 
components. This is so because the gauge fields can , and do, couple the different fermionic 
components. In this fashion we will be left with only slowly varying fields. Thus, with 
the same notation used for the fermions, out of the components of the gauge field A^ we 
define four sublattice amphtudes ^°(x, t) (a = 1, ... ,4). Since the space components Aj 
only enter in the action in exponential hopping amplitudes, it is convenient to define the 
sublattice amplitudes ^^^"(x, t) = exp(i^^(x, t)). In the uniform flux phase {i.e. for A < Ac) 
the hopping amplitudes take the expectation values W"" (j = 1,2). The AFA equations 
tell us that, in the flux phase, the oriented product of the amplitudes W"- around each 
plaquette must be equal to —1. In section |V|we solved this requirement with the choice of 
eq. (|4.15| ). For the purposes of taking the continuum limit, we will choose instead W^ = i 
and W2 = i{—l)°'. The two configurations are related by a gauge transformation and are 
equivalent modulo a flux of 2tt. The time components have zero average. It is natural 
to define a set of slow, fluctuating, fields A^ by identifying W°'{'x,t) -^ W^ exp{iA'^{'x,t)) 
(notice that the A'j is now a fluctuation!). 

These fluctuating fields are small but not slow, that is the fluctuations involve low- 
frequency processes with wavevectors which are not necessarily small. In fact, the gauge 
fields have components with wave vectors which mix all the fermionic components. For this 
reason, we proceed to define a set of fields which are slow and for which there is a sensible 
long distance limit. This we do for both fermions and gauge fields. The sublattice fermion 
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amplitudes can be combined into two species of Dirac spinors \E'^ (labelled by r = 1, 2 , each 
with two components labelled by a = 1, 2), defined by the following linear combinations 

^1 = 7S— (-#1 + #2 + ^3 - ^4) 
Zflo 

*2 = 7^(^l-^2-#3 + #4) (5.1) 

ZOo 

where Oq is a lattice spacing. The normalization factor is chosen so that the term in the 
action which includes the time derivative has coefficient one in the continuum limit. The 
Dirac structure is a consequence of the spectrum of the fiux phase which is linear. Unlike 
the amplitudes ipa which are dimensionless, the continuum fermi fields \1/^ have dimensions 
of (length) ^^. This is the correct canonical dimension for a Dirac field in 2 + 1 dimensions. 
We also need to define a set of gamma matrices which in 2 + 1 dimensions are 2x2 matrices 
which act on the Dirac components of the spinors, labelled by the greek index a. We choose 
them to be 70 = 0"2, 71 = icri and 72 = icr^, where ctj (? = 1, 2, 3) are the 2x2 Pauli matrices. 
We will also need a second set of Pauli matrices, that we denote by T^ (6 = 1, 2, 3), which 
will act on the species index of the spinors, labelled by the latin index r. When necessary, we 
will also use the identity matrix / for each set of indices. In order to simplify the notation, 
we will avoid the explicit use of the indices and, instead, use the matrices to denote the 
operators of interest. 

With these definitions, and after rescaling the time coordinate t by the Fermi velocity vp 
as t = Xq/vf (in units of J, we get v-p = ao! here v-p is measured in units in which the XY 
term of the Hamiltonian has amplitude |), the free part of the fermion action becomes, in 
the continuum limit 

5^°) = f d^X^ri^'^r (5.2) 



with \1/ = ^1^*70. This continuum theory is valid for fiuctuations with wavevectors smaller 
than some cutoff A ^ — , where Oq is the lattice spacing. The fact that this is an effective 
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theory at length scales long compared with the lattice spacing will have important conse- 
quences for our analysis. In field theory the choice of cutoff (or regularization) is largely 
arbitrary. Different choices of regularization usually lead to the same theory. However, for 
theories of fermions a number of subtleties arise connected with the way regularizations treat 
the symmetries of the effective continuum theories. In 2 + 1 dimensions relativistic massive 
fermions have a parity anomaly. However, if the fermion mass has a dynamical origin {i.e. if 
it is induced by fiuctuations) the parity anomaly may be lost in some regularization meth- 
ods (such as dimensional regularization) which set to zero all non logarithmic divergencies. 
In any case, in our problem we are not free to choose an arbitrary regularization scheme 
to cutoff the divergencies present in various Feynman diagrams of the effective continuum 
theory. Instead the choice of cutoff will be done in such a way that the symmetries of the 
lattice system are respected. In particular, we will see in the next section that schemes such 
as dimensional regularization cannot be used for this problem. 

Next we define new linear combinations of the gauge fields. We choose the new linear 
combinations which have a dominant wavevector {qi,q2) such that it mixes the fermionic 
components defined above. We use the notation A'j^'^^. The important wavevectors are 
(^1) Q2) = (0, 0), (tt, 0), (tt, tt), (0, tt). The new fields are 



A, = ^^iAl + Al + Al + Al) (5.3) 

Af = ^^{Al-Al + Al-Al) (5.4) 

A;- = ^^{Al-Al-Al + Al) (5.5) 

Af = ^^{Al + Al-Al-Al) (5.6) 

where A^ = A°° are the smooth fiuctuations of the gauge field. 

The only fiuctuations that are usually kept in these type of analysis are the smooth fields 
like A^. However, we find that some of the other amplitudes are very important. In fact 
we will find that the fiuctuations AZ"^ are connected with Neel fluctuations. The fields Af^'^'^ 
couple to operators which are bilinear in Fermi fields and whose characteristic wavevector 
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is (gi,g2). 

In terms of the fermion amplitudes defined in eq. ( [5 .11) there is a total of 16 operators 
which are bilinear combinations of the amplitudes \1/J^. They have the form M = 'J'\l/, 
M'^ = ^T''\1/, Mfj_ = ^7^\1' and M^ = ^7^T''\E'. In spin language these fermion bilinears 
correspond to linear combinations of site occupation numbers and of hopping amplitudes 
among sites on the four sublattices. Out of the 16 bilinears we will only discuss three of 
them, M, M^ and Mq, which correspond to the parity breaking mass operator, the Neel order 
parameter and the charge density operator respectively. We find the following identifications 

M = ^^^ ^{^Pl^Pi - ^pl^Pl) + ^(^|V'2 - V-I^s) (5.7) 

Mo = ^70* ^ ^Ui + ^y2 + V'sV's + V'4^4 (5.8) 

M^ = ^T^^iJ ^ -^1^1 + ^t^2 + ^1^3 - ^t^4 (5.9) 

These identifications show that Mq is the total fermion occupation number averaged over 
the four sublattices. In the sector with 5*2 = we expect to find (Mq) = 0. The operator M^ 
is the staggered occupation fermion number which, back in spin language is the Neel order 
parameter. Finally, M is an operator which induces hopping across the main diagonals of 
the plaquettes. It is also easy to show that the phase factors present in the definition of 
M, when combined with the phase factors of the fiux phase, indicate that the fiux on every 
elementary triangle inscribed in each plaquette is equal to |. Thus, a non-zero value of (M) 
in the ground state breaks both T and P. M is the Chiral Order parameter introduced 
by Wen, Wilczek and Zee |^3[. We will see in section|V| that this operator plays a very 



important role in the determination of the physically correct (PC) phase diagram. 

The next step is to write the action of eq. ( p.lOp in terms of the slowly varying fields \I^^ 
and AJ^"^^. We discuss first the Fermion part of the effective Lagrangian Cp in the continuum 
limit. After expanding the hopping amplitudes up to leading order in fluctuations and after 
taking the continuum limit we get 

Cf = ^ ip "^ + A^^'^T^'i' (5.10) 
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where p = D^'-/'^ and D^ = d^~iA^ is the covariant derivative which represents the couphng 
to the long wavelength smooth pieces of the Chern-Simons gauge field. We have not included 
in the final form of the Lagrangian, additional terms of the form £extra, 

Axtra = A°-^T27i^ + Af^T^^^^ - Al^^<T^^^^< - Af ^TSo^ (5.11) 

which couple the fermions to fiuctuations with wavevectors (0, vr) and (vr, 0), since it is 
possible to show that they are irrelevant in the renormalization group sense. Physically, this 
can be understood since these operators do not acquire an expectation value in any of the 
phases of our interest. Using the linear combinations of eq. ( ^.1|) , one can show that these 
operators are related with both spin density wave order parameters with wave vectors (Ovr) 
and (vrO) and to Peierls (or dimer) order parameters with the same wave vectors. Neither type 
of order occurs in this system. From now on we will ignore these terms. Notice that these 
terms have the same scaling dimension as the ones that are being kept. It is conceivable that 
there are other situations in which these terms become dominant such as in the vicinity of a 
dimerization transition. Such phase transitions are possible in a frustrated antiferromagnet. 
However, the Chern-Simons ( or Wigner- Jordan) mapping used here cannot be used as it 
stands for non-bipartite lattices. The terms in £extra break the rotational invariance of the 
continuum theory of eq. (|5.10|) down to the allowed symmetries of the point group of the 



square lattice. Again, such anisotropies will only be relevant in the vicinity of ground states 
which break such symmetries. Below we will also ignore terms in the bosonic part of the 
Lagrangian with the same symmetry properties. 

Finally, we need to find the continuum form of the bosonic parts of the action ( p. 101) . 



We only keep terms to leading order in the lattice spacing for each field. There are two sets 
of contributions. One set comes from Si^f When written in terms of the slow components 
of the gauge fields they contribute with the term £int 

Ant = I {AY - AYf - ^ {d,A2 - d2A,f (5.12) 

where the effective coupling constants g and e^ are 
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4^ ^2 _ ^2 (^g_^3) 



Q = ^ — e 



Xao 



The Maxwell-like term in the Lagrangian Ci^t has no consequence on the phase diagram 
since it has one more derivative than the Chern-Simons term and, hence, it is irrelevant at 
long distances |^ and it will be ignored from now on. 



The second set of contributions comes from the Chern-Simons terms Scs- We find a 
contribution to the Lagrangian of the form 

n no 

£cs = Te'^"" A^F.A + - AJ" {AT - AY) (5. 14) 

4 ao 

where F^^ = dfj,Ai, — d^A^ is the field strength. 

Before putting everything together we notice, by inspecting eqs. ( |5 . 1 2|) and (|5.14|) , that 



the fields A^ and A^ enter in the total Lagrangian in terms of the difference A^ — A^ 
and that the total Lagrangian is at most quadratic in this difference. Thus, we can integrate 
out these fields and find the Lagrangian of the form 

£ = ^ ^ + AY^T''^ - ^ {A^f + -f^^'^A^F.x (5.15) 

Zg 4 

where g is an effective coupling constant. It is convenient to define a dimensionless coupling 
constant (70 such that g = gocto- In terms of the anisotropy A and of the Chern-Simons 
coupling 6 we get 

where we have introduced the statistical angle 6 = 1/26. Notice that for ^ ^ ^^ {i.e. in the 
boson limit of interest here) 5 — > tt and go —>■ 4X. Conversely, in the fermion limit 6 —>■ or 
6 ^ 00, the dimensionless coupling constant is weak, g^ —^ 0. Thus, in the hard-core boson 
limit, which is the case of interest, we are dealing with a system in which the dimensionless 
coupling constant is typically of order unity and perturbation theory should not be reliable. 
Let us discuss the physical meaning of the terms of the Lagrangian C of eq. (p.lSp. 



The field A^ couples locally to the Neel order parameter operator M'^ = ^T^^. Thus, if 
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Aq'^ picks up an expectation value, so does the order parameter M^. Clearly such a state 
has Neel long range order. In section |IV| we found that beyond some critical value of the 
anisotropy A the system is in a Neel state. The field Aq'^ also enters at most quadratically in 
the Lagrangian C and it can also be integrated out (in fact, we may regard the field A'^'^ as 
a Hubbard- Stratonovich field). Indeed, after integrating out Aq^ the total Lagrangian £20 
for the two-dimensional system has the suggestive form 

£2D = ^ ^ + I {^T^^) ' + ^e^^'^A^F^x (5. 17) 

Thus, we find that the low energy degrees of freedom of the anisotropic quantum antifer- 
romagnet can be described in terms of a theory of two relativistic Fermi fields in 2 + 1 
dimensions, with a four-Fermi interaction of strength g, and which are also coupled with a 
Chern-Simons gauge field. The Chern-Simons coupling constant is restricted to the value 

= ^. For the rest of this paper, we will consider this Lagrangian. 

In the form of eq. ( p.l7| ), the effective continuum theory for the two dimensional system is 
a generalization of the equivalence that exists between the one-dimensional anisotropic spin- 

1 quantum Heisenberg antiferromagnet and a theory of a self-interacting relativistic Fermi 
field, the Luttinger-Thirring model PJIDj] . The left and right components of the Fermi fields. 



^E^L and \l/pi,, can be viewed as the two components of a Dirac spinor in 1 -|- 1 dimensions. 



The Lagrangian Cid of the effective field theory for the antiferromagnetic spin chain is |44 



£iD = ^i^^ + — (^^)^ (5.18) 

The coupling constant for the ID system is qid = 2A which is dimensionless. 

The two-dimensional theory, with Lagrangian C21), and the one-dimensional theory with 
Lagrangian Cid differ in a number of ways: (a) in ID there is only one species of Dirac 
fermions (instead of the two labelled by r in £20 ), (b) there is no gauge field in ID but 
there is a Chern-Simons gauge field in 2D, (c) cjid is dimensionless while g has dimensions 
of length and (d) the symmetries are different. They both have a self-interacting, quartic, 
term in fermions which, in both cases, is equal to the square of the Neel order parameter. 

32 



We will also see that the coupling to the Chern-Simons term is, in this problem, responsible 
for much more than a change in statistics. 

VI. CRITICAL BEHAVIOR AND PHASE DIAGRAM 

In this section we discuss the phase diagram of the effective field theory derived in 
section^ with the Lagrangian £20 of Eq. (|5.17|) . Given the apparent similarity between 



this theory and its analog in one dimension £id of Eq. ( ^.1^ ), one might think that the 
phase diagrams may also be quite similar. However, a closer analysis shows that this cannot 
possibly be correct. 

The two theories have different symmetries as well as different scaling properties. In 
both cases, the order parameter of the Neel state is odd under a sublattice exchange or, 
what is the same, it is odd under a global shift of the field configuration by one lattice 
spacing. This is a global discrete symmetry which can be spontaneously broken by the 
ground state, even in one space dimension. In the case of the one dimensional chain, the 
Neel order parameter, i.e. the difference of the spinless fermion occupancy between the two 
sublattices, is proportional to ^\1/. This operator is odd under the transformation \1/ -^ 75 \l/ 
which changes the relative sign of the right and left moving amplitudes of the fermions. This 
symmetry is known as a discrete chiral symmetry. The Lagrangian Cid is even under this 
chiral symmetry. As a result, it is possible to show that the operator ^\I^ does not acquire 
an expectation value to all orders in perturbation theory in the coupling constant guj. Also, 
because this symmetry is present, renormalization effects do not induce fermion mass terms 
in the Lagrangian, which are proportional to 'l'\E' and thus brake the symmetry explicitly. 

The 1 + 1-dimensional system has very special scaling properties. In space-time dimen- 
sions D = 2, the coupling constant qid is dimensionless. Then, the standard field-theoretic 
analysis tells us that the theory defined by the Lagrangian £id, the Gross- Neveu model |^ 



is renormalizable. If the number of fermionic species is A^, this theory is asymptotically free 
with a /3-function [Q Pigin) = i^^) 9id- Fo^ N > 2 this /3-function is strictly positive 
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and the system has a dynamically generated energy gap. For the case of interest for the one- 
dimensional quantum Heisenberg antiferromagnet, A^ = 1 and the coupling constant is still 
dimensionless. However, in this case, the /5-function vanishes to all orders in perturbation 
theory and the system has a line of fixed points. In other words, the four fermion operator 
(i|rv[/)2 jg marginal. Bosonization studies show P, ^^ | that there is an operator which is 



irrelevant at small coupling but that at a critical value of the coupling constant (namely 
for a critical anisotropy) it becomes relevant. The effect of this operator, which represents 
umklapp processes in the lattice theory, is to end the line of fixed points at a multicritical 
point which is in the Kosterlitz-Thouless universality class. An explicit computation of the 
critical exponents shows that the correlation functions at this multicritical point have 
the symmetries of the isotropic Heisenberg antiferromagnet. 

In the case of the two-dimensional system, the order parameter is 'l'T^\I'. Unlike one 
dimension, there is no 75 Dirac matrix in D = 2 + 1. Instead, the order parameter is now 
odd under an operator which effectively exchanges the two species of fermions. However, a 
mass term proportional to the operator 'l'\E' is even under the exchange of fermionic species 
(or sublattices, which is the same) but it is odd under parity (P) and time reversal (T) 
transformations. In contrast, ^T^'^ does not break these symmetries. The reason is that, 
for each species, the operator ^\1/ changes sign under P and T and, hence, ^T^'^ changes 
sign too. However, this effect is equivalent to a redefinition of the sublattices and, thus, it 
breaks neither parity nor time reversal. For this reason, the operator i|'T^\l/ is usually called 
a parity invariant mass term. Thus, renormalization effects cannot induce a parity invariant 
mass term in perturbation theory since it is a symmetry breaking field. We will show below 
that the phase transition found in section |I^ represents the spontaneous breaking of the 
symmetry ^T^"^ — > -^T'^^ to a Neel phase in which (^T^\&) ^ 0. 

Unlike the 1 + 1-dimensional theory, the Lagrangian £20 of eq. ( |5.17| ) has, in addition 



to the discrete symmetry ^T^^ -^ \I^T^\I', a continuous gauge (local) symmetry, 

*(x) -^ e^'^(^)^(x) 
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A^{x) ^ A^{x) + 9^0(x) (6.1) 

With the rules that we are using here, the gauge field has scaling dimension one and the 
Chern-Simons term has scaling dimension three and it is marginal. The Chern-Simons 
coupling constant 6 is dimensionless and the theory is renormalizable. In contrast, the four 
fermion operator has naive scaling dimension four and it is irrelevant. The four fermion 
coupling constant g has naive scaling dimension of length. The dimensionless coupling 
constant ^^o of eq. (|5.16|) has a negative /3-function, P{go) = —go + O^g^) and (^g = is 
an infrared stable fixed point of the renormalization group fiow. In contrast, the Chern- 
Simons coupling 6 has a vanishing /3-function. Thus, at least naively, it appears that this 
renormalization group fiow has a line of fixed points parametrized by 6. 



W. Chen and M. Li ^^ have recently put forward arguments in favor of a scenario in 
which theories of relativistic fermions coupled to Chern-Simons gauge fields have a line of 
fixed points parametrized by 9. In particular they argue that the fiuctuations of the gauge 
fields induce anomalous dimensions in the fermion operators and, as a result, the four- 
fermi operators become less irrelevant as 6 grows larger. However, their analysis is based in 
dimensional regularization. It is well known that this regularization scheme sets to zero all 
singular Feynman diagrams except for logarithmically divergent terms. Since it is an analytic 
regularization method, it also sets to zero all contributions which cannot be continued in 
dimension. Thus, fermion mass terms are not induced in dimensional regularization. 

At the infrared stable fixed point g = 0, the Lagrangian £20 is manifestly scale invariant. 
The only parameter left is the Chern-Simons coupling constant 6 which is dimensionless. 
Hence, at least in the absence of fiuctuations, it is at a fixed point and it has no scale. 
However, it is easy to convince oneself that fiuctuations lead to divergent corrections of 
the classical (that is, mean field) values of the observables. The presence of divergent 
contributions in the perturbation series requires the use of a cutoff or, more generally, of a 
regulator. Any cutoff introduces a microscopic scale in the problem and therefore it breaks 
the apparent scale invariance. Thus, in the presence of a cutoff, dimensionful terms can 
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be induced by renormalization effects. On dimensional grounds, an induced mass term will 
have to be proportional to the cutoff since there is no other scale left. The theory that we 
are studying has a natural cutoff, the lattice spacing, which acts as a natural scale. Hence, 
in this case, it is physically incorrect to use dimensional regularization. Furthermore, the 
Chern-Simons term breaks both T and P. Thus, it is expected that renormalization effects 
should generate all possible terms which break the same symmetries. The parity-odd fermion 
mass term ^\1/ breaks precisely the same symmetries. Hence, unlike the 1 + 1-dimensional 
case there is no symmetry that prohibits these terms to be induced by renormalization. 

Unlike the 1 + 1-dimensional theory, non-perturbative tools such as bosonization are not 
available for the study of relativistic systems in 2 + 1 dimensions. In order to proceed further 



we will use a semiclassical approximation, like the one of section |V|, but going beyond 
the leading order. In 1 -|- 1-dimensions this approach would not be sufficient since this 
approximation misses the marginality of the interaction. However, in 2 + 1 dimensions the 
four fermion interaction is irrelevant (in weak coupling) and this approximation reproduces 
this result correctly. In order to determine the induced fermion mass Mind, we will compute 
the leading self-energy correction to the fermion propagator due to fluctuations of the gauge 
field. Since the gauge field couples in the same way to both species of fermions, the only 
possible induced mass term is ^'^ which is even under the exchange of species. Since this 
term is odd under T and P it can only arise from fluctuations of the Chern-Simons gauge 
field. 

The partition function at zero temperature is 

where C is the Lagrangian of Eq. (|5.15| ), which we reproduce here for clarity 

1 f) 

C = ^ ip-^ + Al'^^lT^m - — {Al'^f + -e'^^M^F^A (6.3) 

Zg 4 

This form of the theory makes the semiclassical approximation more transparent. We now 
follow the methods outlined in section ^ and find an effective action S'eg for the Bose fields, 
which here are A^ and Aq'^ , after integrating out the fermions: 
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ip + Al^T'] -Jd'x— (A^-)' + J d'x -e^'-'A.F.x (6.4) 
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The stationary points of this action satisfy the equations 



^^ ^^^ t Tr [S{x, xh,] + -e^.xF'-' = (6.5) 



5S, 



eS 



i Tr 



S{x,x)T^]--A^''{x) = (6.6) 



M5^(x 
where S'(x, x') is the (Feynman) fermion propagator 

S{x,x') = {x\- — ^ :r^^W) (6-7) 

In order to simphfy the notation we have dropped all the indices that are attached to the 
fermions. The traces run over both spinor (Dirac) and "flavor" (species) indices. 

The equations ( |6.6|) have a solution of the form A^ = and A^^ = M. Please recall that 
the gauge field A^ is the long wavelength fluctuation of the (lattice) Chern-Simons gauge 
field around the flux phase. M is the Neel mass and it is given by 



g J (27r)3 ^-MT^ + ie' 



which is the continuum analog of the gap equation of section |I^. After performing the 



integral (with a cutoff on the space components of the momentum A ^ 2^) we find that the 
mass M is the solution of 



M A / M2 IMI \ , , 

This gap equation has the solution 

m = U^-'-]Qi9-9c) (6.10) 

^ \9c 9 J 

where B(t) is the Heavyside function and cjc = j^ is the critical coupling constant of the 
effective continuum theory. Using this value of cjc and eq. ( [5.16| ) we get an estimate for 
the critical anisotropy Ac from this continuum theory given approximately by Ac ~ ^ which 
should be compared with the (lattice) AFA value Ac ~ 0.4 of section |^. This different 

37 



value of Ac, which is a non-universal quantity, reflects the approximations made in taking 
the continuum limit. In particular, it depends on the precise relation between the momentum 
cutoff A and the lattice spacing a^. 

Let us denote the full fermion propagator at momentum p^ {fi = 0, 1, 2) by S{p) and by 
G^'^{p) the propagator of the Chern-Simons gauge field. At this level of approximation, the 
fermion propagator, at 3- momentum p^, is given by 

Soip) = ; \,, ^ . (6.11) 

where, once again, we have dropped all indices. The bare propagator of the gauge field at 
3-momentum p^, in the Lorentz gauge d^A^ = 0, is 

p^ +ze 
Thus, in every Feynman diagram, each propagator of the gauge field contributes with a 
weight proportional to 1/6 = 26, where 6 is the statistical angle. Therefore, this is an 
expansion in powers of the statistics and it is accurate only near the fermion limit 6 ^ or 
6 ^ oo. The value of 6 of interest for the Heisenberg antiferromagnet is 5 = 27r which is not 
small. A more serious problem is that the physical properties of systems of this sort must 
be periodic in the statistics. Namely, all amplitudes for any physical observables must not 
change under the replacements 6^6 + 2Tfk, where k is an even integer (periodicity) and 
6 -^ 4t[ — 6 (symmetry around bosons). We will use parturbation theory around fermions 
and demand that it holds around each period. The extrapolation to the boson point 6 = 27r 
should yield qualitatively correct results. 

The exact fermion propagator S{p^) obeys the Dyson equation 

S{p)-' = So{p)-'-^{p) (6.13) 

where S(p) is the fermion self-energy. Due to the symmetry of the bare theory, S(p) has 
the form 

J:rr'ip) = ^lip) + MP)T' (6.14) 
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To leading order in S, S(p) is given by 



S(p) 



Explicitly, we find 
S(p) = 



d^k 



(2vr) = 



Y 



rSoip-k)r^G^o''ik) +0(6' 



(6.15) 



^-]^ + MT^ 



-1 -^(^tJ.u\- 



le 



0{5^ 



(6.16) 



(27r)3 ' (p-A;)2-M2 + ze ' 6 ^'^^ k^ + ie 
By counting powers of the momentum of integration we see that this contribution has an 
ultraviolet, linear divergence. By expanding the self energy in powers of the external mo- 
mentum p^, we see that the (ultraviolet) linear divergence only affects the term at zero 
external momentum p = and that all contributions at non-zero external momentum are 
finite. This ultraviolet divergence is an artifact of ignoring the fact that we are working with 
an effective field theory and that the quantum antiferromagnet, from which this field theory 
is derived, is defined on a lattice and it does have a cutoff. Thus, this divergence has to be 
cutoff at values of the internal momentum of the order of A ^ 7r/2ao , where oq is the lattice 
spacing. After some algebra we find 



E(p) 



^dx ( "^'^ {qx + xpx)[{p,{l~x)~qp)2g{-2MT^^y\ 



(6.17) 



e Jo '^"^ J (27r)3 [g2 + ^23,(^1 - s) - M'^xY 

where we have used the covariant notation q^ = q^ — q^ and an ie prescription is assumed. 
After an integration over the frequency variable go and over the spacial components of 
the momentum of integration q (with a cutoff A) we get 



s(rt 



;dx^(xMn,,^, 



1 



1 



M'^x — p^xil — x) \ K^ + M^x — p^xil — x) 



2n9 

1 
"4^ 



-i/A^ + M'^x — p'^x{l — x) — J M'^x — p^x{\ — x) 



f— - (m^x - 2x(l - x)/) 



Ym^x - p%(l - x) ^ K^ + M^x - p'^x{l - x) 
Since we are only interested in the computation of the effective (or renormalized) mass, it 
will be sufficient for our purposes to compute the integral of eq. ( |6.18| ) at p = 0. In this limit 
we find that S2(0) = and we get a value for S(0) which is independent of the fermionic 
species. In this limit, and after some algebra, eq. (|6.18|) becomes 



(6.18) 
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This contribution to the self energy of the fermion plays the role of an effective or induced 
mass and we will denote it by S(0) = Mjnd- Due to the symmetry of the bare fermion 
propagator, this (divergent) induced mass is the same for all fermionic species which thereby 
acquire the same mass. This result also tells us that the induced mass is proportional to 
— 1/6. The fact that this sign is opposite to the sign 6 will play a fundamental role in our 
analysis. 

Let us now use these results to compute the mass of each species of fermions up to 
corrections of order 6"^. Our calculation tells us that the total fermion propagator S{p^) at 
zero external momentum has the form 

S-i(O) = -MT3-S(0) (6.20) 

from where we find that the effective masses Mi {i = 1, 2) for each of the species are 



--<-)--^U-f -Tl 



where M is the solution of the gap equation ( |0| ) and it is a function of the coupling constant 
g. Hence, the effective masses Mj are also functions of the coupling constant. In Fig. [I| we 
show the qualitative form of the functions Mi{g) for the entire range of couplings. Given 
the relation between g and the anisotropy A, the curves of Fig. || will help us to determine 
the phase diagram. 

Thus, while the spectrum of the semiclassical theory consists of two massless fermions 
for all g < cjc and two massive fermions (but whose masses have opposite signs) for g > cjc, 
the quantum fluctuations of the gauge fields make the fermion spectrum generically massive 
for a// values oi g ( and, hence for all A). However, eq. ( |6.21| ) shows that, for 6 positive and 
9 < dc, Ml = M2 = —^ < 0. In fact, for 6 > 0, M2 is always negative. However, Mi goes 
through zero and changes sign at some critical value of the coupling constant g*{0) 



r{e) = -gJl + j^ (6.22) 
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For the value oi 9 = l/27r, of interest for the Heisenberg antiferro magnet, we get g* = -\/3^c- 
Thus, we conclude that the quantum fluctuations yield the following spectrum for the 
fermions. For g < g* the two species of fermions have masses which, for general values of 
g have different absolute values but have the same sign. This sign is opposite to the sign 
of 6. In contrast, for g > g* the two fermionic species have masses with different absolute 
values and opposite signs. Precisely at g*, the mass of one of the species of fermions passes 
through zero. We should regard this phenomenon as a phase transition and g* as a critical 
point. From the arguments presented above we should expect that the critical value g* 
should correspond to a critical value of A, which we will denote by A*, and that there should 



be a phase transition at A* in the anisotropic quantum Heisenberg antiferromagnet ^1 

Given that the fluctuations make such important effects already at the level of the leading 
corrections to the AFA, it is natural to inquire what are the effects of even higher order 
corrections. It is clear that nothing special happens at cjc and that it does not correspond 
to a phase transition which has been shifted to g*. The apparent discontinuity in the 
derivative of the masses at cjc is an artifact of the leading order calculation. Higher order 
terms will smooth out this spurious effect. The actual value of the critical coupling will also 
be renormalized by higher order terms. 

We now will argue that the phase transition a.t g = g* should be identified with the 
isotropic Heisenberg antiferromagnet. The argument in support of this identification relies 
on the counting of massless excitations of the system. In turn, we will identify these massless 
excitations with the Goldstone modes of the antiferromagnet. 

At this level of approximation, the spectrum of fermions is massive for generic values 
of the anisotropy with one of the species becoming massless just at the critical coupling. 
Let us now investigate the spectrum of the bosonic excitations, A^ and A^'^, as a function 
of g. Since the fermions are massive, it is possible to integrate them out off the partition 
function and to find an effective action for the Bose fields which is local at length scales long 
compared with the inverse of the mass. From the work of Deser, Jackiw and Templeton ||2^ 
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we know that the long distance effective Lagrangian for the gauge field, £ind(^^), induced 
by the fluctuations of a fermion of mass M, is of the form 

C-UA,) ^ -^F^-'F,, + ^-^^^e.^.A^F'"' (6.23) 

where we recognize the last term as a Chern-Simons term and the sign of its coupling 
constant is equal to the sign of the fermion mass M. The parameter k is proportional to 
\M\. A similar analysis implies that the fluctuations of A^'^ are always massive. However, 
since the fermions are actually massive at Qc, the fluctuations of the Neel order parameter 
field Aq"^ never become critical. 

The total effective action for the fluctuating gauge field is the sum of three contributions: 
(a) the (bare) Chern-Simons term with coupling constant 6, and (b) two additional Chern- 
Simons terms each with the value of eq. (|6.23|) . Thus, the total Chern-Simons coupling 
constant, O^s is 

9,,te).9H.!!i£WM±i§M (6.24) 

4:TC 



eos{-g) = { '^ (6.25) 



We find, 

' e-l \i-g<r 

6 otherwise 

Thus, in the phase in which the effective fermion masses have the same sign, the fluctuations 
of the fermions act so as to reduce (or screen) the bare value of the Chern-Simons coupling 
constant 9. This is the "Lenz law of statistics" referred to in the Introduction. 

For the particular case of interest for the antiferromagnet, 9 = 1/2tt, the screening is 
complete and ^cff = for the case of bosons!. Thus, the Chern-Simons term is cancelled out 
from the effective action of the gauge field which now represents massless excitations. This 
cancelation of the Chern-Simons term from the effective action of the gauge field is the anyon 



superfluid scenario of references p0|-|2^. In the anyon superfluid, the transverse massless 
gauge field is interpreted as the Goldstone boson of the superfluid state. We identify this 
regime with the XY phase of the antiferromagnet. 
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Conversely, for g > g*, we get ^eflf = & and all collective modes are massive. In this phase 
the Neel order parameter is non zero and all excitations have a gap. This is the Ising phase 
of the antiferromagnet. 

The phase diagram is now almost complete. Two issues still need to be resolved: (1) the 
nature of the transition point at g* and (2) are the fermion states really part of the spectrum 
for any value of the coupling constant? 

These two problems actually are not independent from each other. Let us first consider 
the fate of the fermions. For g > g*, the fermions have non- zero masses with opposite 
sign. In this phase, the effective Chern-Simons coupling constant is equal to l/27r. Hence, 
in this phase, by the standard argument of statistical transmutation, Chern-Simons gauge 
field turns the massive fermions into massive bosons. However, for g < g*, the long distance 
effective action for the gauge field does not include a Chern-Simons, which cancels out 
but, instead, the leading term has a Maxwell form. Thus, the actual physical mass (or 
gap) of the fermion will be significantly renormalized by the fluctuations of the transverse 
massless collective mode (the gauge field). The result, exactly as in the case of the semion 
superfluid [^ , is that due to the quantum fluctuations of the collective mode, the fermion 
acquires a logarithmically divergent mass and, therefore, it disappears from the physical 
spectrum. 

The presence of infrared divergent corrections to the fermion self-energy in the range 
g <g* has important consequences for the Neel order. Eq. (|6.6|) relates the expectation value 
of the Neel order parameter '|'T^\I' to the expectation value of the field Aq'" (this equation is 
valid beyond the saddle point approximation provided that the field Aq'" is replaced by its 
exact expectation value). At the level of the saddle-point, the expectation value of Aq'^ is 
equal to Mg. Higher order corrections involve fermion self-energy insertions in the low order 
diagrams. For g > g*, these corrections are finite in the infrared. However, for g < g*, these 
corrections are infrared divergent (due to the fermion self-energy insertions) and negative 
(since they have to give a value smaller than the classical result). The precise computation 
of these effects to all orders is difficult and beyond the domain of perturbation theory. 
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The presence of these infrared divergent contributions, aheady in the leading corrections, 
suggests that the Neel order is unstable for g < g* and that the exact expectation value 
of this operator has to vanish in that regime. An alternative picture of this effect can be 
seen by noticing that, as the coupling constant g is decreased from large values ( that is 
from the classical Neel regime) and one of the masses becomes small, tunneling processes 
between the two Neel states become increasingly favorable. In particular, the magnitude 
of the energy per unit length of a Neel domain wall is set by the mass of the fermionic 
excitations. Thus, in the regime in which the fermionic excitations have a finite mass, the 
energy per unit length of the domain wall is finite. Below g*, the infrared divergencies in the 
fermion self-energy will force the energy per unit length of the wall to vanish. Consequently, 
the domain walls will condense in this regime and will destroy the Neel long range order. 
Precisely at g*, the energy per unit length of the domain wall is still finite since there is 
still Neel order. Thus, the domain wall energy should drop to zero with a jump at g*. The 
domain wall condensation as a mechanism for the destruction of long range order is well 
known in 1 + 1-dimensional systems, where the solitons play the role of the domain wall. 
Hence, we conclude that the Neel order parameter should drop to zero discontinuously at 
g* and to vanish for all g < g*. 

In contrast, the operator which creates (or removing) one fermion and one flux quantum 
simultaneously, is gauge invariant and it has a finite mass. In anyon superfluidity this state 
is usually called the vortex and it has the statistics of an anyon. In our problem, again by 
statistical transmutation, it is a massive boson. The mass of this state should scale with Mi. 
Hence, it should vanish exactly at g*. In the lattice Chern-Simons theory this state is created 
by an operator which changes both charge and flux. For the case of the antiferro magnet, 
the operator which creates this state is S~^. Thus, we argue that the massless fermion of 
our spectrum is actually the extra Goldstone boson of the antiferro magnet. Since this state 
becomes massless only at g*, we identify this phase transition with the isotropic quantum 
antiferromagnet. This point is then viewed as a limiting point and it has the attributes of 
both phases. In particular, it has a non-vanishing value for the Neel order parameter and 
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two Goldstone bosons 0]. We also note, in passing, that in one- dimension a strikingly 



similar picture of the nature of the isotropic point (from the point of view of the symmetry 
analysis) was developed by Luther and Peschel 0. 

However, in the case of our problem, unlike the case of the one-dimensional spin chain, 
there are no non-perturbative methods available, such as bosonization, that will allow for 
an exact treatment of the long distance behavior. The topological invariance of the Chern- 
Simons theory strongly suggests that no further renormalizations occur and that we have 
successfully characterized the infrared stable fixed points. A detailed analysis of the phase 
transition at g*, particularly the determination of its universality class, critical exponents, 
etc., requires a more sophisticated analysis of the effective theory than the one performed 
here. The phase transition at cjc is in the universality class of theories of self-interacting 
relativistic fermions (of the Gross-Neveu type). Our analysis shows that this fixed point is 
unstable and does not represent the long distance behavior of the system. It is interesting 
to note that the removal of the phase transition at cjc and its replacement by the transition 
at g* where the Neel order parameter drops to zero dis continuously, is strongly reminiscent 
of the physics of fluctuation induced first order transitions. 

We can also give a renormalization group picture for the arguments presented above. The 
AFA or, equivalently, the semiclassical theory of this section, does not represent faithfully 
all the relevant fluctuations of the system. In particular, the infrared stable fixed point 
associated with the XY phase is simply not present in the semiclassical theory. The quantum 
fluctuations of the gauge field contain the appropriate relevant operator, the parity breaking 
induced mass term. Once this operator is generated, the flow of the coupling constants is 
drastically changed as we explore the low energy regime. In particular, the effective action 
of the gauge fields acquires a finite renormalization of its Chern-Simons coupling constant 
which tends to screen the statistics. The fermions are also affected by this flow since the 
physical low energy mass (or energy gap) of the spectrum is altered by the fluctuations of 
the bosons. In the XY phase, the excitations of the gauge fields are massless and their 
quantum fluctuations suppress the fermions from the spectrum. In the Ising phase, their 
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fluctuations turn the fermions into bosons. Clearly at both stable fixed points, parity and 
time reversal are not broken (for 6 = l/27r). 

We conclude this section with a few remarks. Firstly, our arguments show that, for all 
values of the anisotropy, there are no states in the spectrum with the quantum numbers of a 
fermion. All the states are bosonic. This is not an accident since the system is not in a spin 
liquid state for all values of the anisotropy. Nevertheless, given the analogies between the 
AFA and the fiux phases of the theories of frustrated antiferromagnets, our results should 
be viewed as an indication that, once fiuctuations are fully taken into account, the fiux 
phases could become more like the standard phases of antiferromagnets. We have presented 
qualitative arguments which show that the phase transition at g* can be viewed, in some 
sense, as a first order transition since the order parameter must have a jump at that point. 
A more rigurous proof of this statement still needs to be constructed. Nevertheless, the 
arguments presented above show clearly the physical mechanism behind this phenomenon. 
A direct computation of the spin correlation function {S^{x)S'{x')) should demonstrate 
quite explicitly the presence of an additional massless state at the critical point. Work on 
this problem is currently in progress. 

VII. CONCLUSIONS 

In this paper we have presented a field theory for the anisotropic quantum antiferro- 
magnet on a square lattice, based on the Chern-Simons or Wigner-Jordan approach. We 
discussed in detail the phase diagram, as a function of anisotropy, at the level of the Average 
Field Approximation. We found a continuous, second order phase transition at a critical 
anisotropy from a fiux-like phase to an Ising phase. We showed that this phase transition 
is spurious and that the massless fiux phase cannot possibly describe the the XY regime of 
the quantum antiferromagnet. We used a semiclassical theory, based in the method of spin 
coherent states, to derive an anisotropic non-linear sigma model for this system and derived 
a phase diagram for the system, valid in the limit of large spin. We considered the role of 
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fluctuations around the AFA and showed that they induce relevant operators, not included 
in the AFA, which drive the low energy behavior of the system. We derived an effective 
field theory of self interacting fermions coupled to Chern-Simons gauge fields and showed 
that its fluctuations contain all the necessary relevant operators to yield a correct phase 
diagram. In particular they induce P and T symmetry breaking fermion mass terms which 
should necessarily be present for arbitrary values of the Chern-Simons coupling constant. 
We gave a set of arguments which indicate that the fluctuations of the gauge fleld drive the 
theory away from the universality class of theories of self-interacting relativistic fermions (of 
the Gross-Neveu type). We identifled the infrared stable flxed points and showed that the 
spectrum of the system at theses flxed points coincides with the expected spectrum of the 
phases of the antiferromagnet discussed in section 0. 

We conclude with a comment on the accuracy of the AFA. The AFA has become a stan- 
dard tool and it is widely used in a variety of flelds of condensed matter, most prominently 
in theories of the Fractional Quantum Hall Effect. The difficulties that we found in applying 
the AFA to the quantum Heisenberg antiferromagnet show that this approximation can be 
unreliable if the resulting fermion spectrum is massless. In such situation fluctuations may 
(and generally will) induce relevant operators which will necessarily generate gaps in the 
fermionic spectrum. This is not a problem for theories of the incompressible states of the 
FQHE but could well be the case for the compressible even denominator states. 

VIII. ACKNOWLEDGEMENTS 

A. G. R. thanks the Department of Physics of the University of Illinois for its kind 
hospitality. This work was done in part while A. G. R. was a member of the Department 
of Physics of the University of Tennessee at Knoxville. This work was supported in part by 
the National Science Foundation through the grants NSF DMR-9015771 at the University of 
Tennessee and NSF STC-9120000 at the University of Chicago (AGR), NSF DMR-91-22385 
at the University of Illinois at Urbana- Champaign (EF), by a Fellowship of the American 



47 



Association of University Women (AL), and by the U.S. Department of Energy through 
contract DE-AC05-84OR21400 administered by Martin Marietta Energy Systems Inc. 



48 



[1] 

[2] 
[3] 
[4] 
[5] 
[6] 
[7 



[9 

[lo: 
[11 

[12 

[13 
[14 

[is: 



REFERENCES 

J. A. Bednorz and K. A. MuUer, Z. Phys. B64, 189 (1986). 

P. W. Anderson, Science 235, 1196 (1987). 

H. Bethe Z. Phys. 71, 205 (1931) 

P. Jordan and E. P. Wigner Z. Phys. 47, 631 (1928). 

E. Lieb, T. Schultz and D. C. Mattis, Ann. Phys. (N. Y. ) 16, 407 (1961). 
T. Holstein and H. Primakoff, Phys. Rev. 58, 1098 (1940). 

C. Kittel, Quantum Theory of Sohds, John Wiley & Sons , 1987 Chapter 4. 

F. D. M. Haldane, Phys. Lett. A93, 464 (1983). 

A. Luther and I. Peschel, Phys. Rev. B12, 3908 (1975). 
F. D. M. Haldane, Phys. Rev. Lett. 50, 1153 (1983). 

E. Manousakis, Rev. Mod. Phys. 63, 1 (1991). 

T. Kennedy, E. H. Lieb and B. S. Shastry, Phys. Rev. Lett. 61, 2582 (1988); K. Kubo, 
and T. Kishi, Phys. Rev. Lett. 61, 2585 (1988). 

F. Dyson, E. H. Lieb and B. Simon, J. Stat. Phys. 18, 335 (1978). 
J. D. Reger and A. P. Young, Phys. Rev. B37, 5024 (1988). 



E. Manousakis, reference 111, and references therein. 



[16] S. Liang, B. Doucot and P. W. Anderson, Phys. Rev. Lett 61, 365 (1988); S. Liang, 
Phys. Rev. Lett. 64, 1597 (1990). 

[17] T. Dombre and N. Read, Phys. Rev. B38, 7181 (1988); E. Fradkin and M. Stone 
Phys. Rev. B38, 7215(1988); F. D. M. Haldane, Phys. Rev. Lett. 61, 1029 (1988); 
L. loffe and A. Larkin, Int. J. Mod. Phys. B2, 203 (1988); X. G. Wen and A. Zee, 

49 



Phys. Rev. Lett. 61, 1025 (1988). 

[18] S. Chakravarty, B. I. Halperin and D. Nelson, Phys. Rev. Lett 60, 1057 (1988). 

[19] E. Fradkin, Phys. Rev. Lett 63, 322 (1989). 

[20] R. B. Laughlin, Phys. Rev. Lett 60, 2677 (1988); C. B. Hanna, R. B. Laughlin and 
A. L. Fetter, Phys. Rev. B40, 8745 (1989). 

[21] Y. -H. Chen, F. Wilczek, E. Witten, and B. I. Halperin, Int. J. Mod. Phys. B B3, 1001 

(1989). 

[22] E. Fradkin, Phys. Rev. B42, 570 (1990). 

[23] S. C. Zhang, T. Hansson and S. Kivelson, Phys. Rev. Lett 62, 82 (1989). 

[24] Ana Lopez and Eduardo Fradkin, Phys. Rev. B44, 5246, (1991). 

[25] R W. Anderson, Science 235, 1196 (1986). 

[26] I. K. Affleck and J. B. Marston, Phys. Rev. B37, 3774 (1988). 

[27] G. B. Kotliar, Phys. Rev. B37, 3664 (1988). 

[28] Y. R. Wang, Phys. Rev. B43, 3786 (1991). 

[29] S. Deser, R. Jackiw and S. Templeton, Phys. Rev. Lett 48, 372 (1982). 

[30] We thank S. Kivelson for suggesting this interpretation. 

[31] This mechanism is somewhat analogous to the metal insulator transition recently studied 
by W. Chen, M. R A. Fisher and Y. S. Wu (CFW), Phys. Rev. B 48 , 13749 (1993). 
Unlike the problem studied by CFW, in the mechanism that we use here the fermionic 
masses are actual fluctuating flelds which, for some range of A, acquire non-vanishing 
expectation values. 

[32] B. I. Halperin, T. Lubensky and S. K. Ma, Phys. Rev. Lett 32, 292 (1974). 



50 



[33 

[34 

[35 
[36 

[37; 

[38 
[39 
[40 

[41 



[42 



[43 



S. Coleman and E. Weinberg, Phys. Rev. B7, 1888 (1973). 

Eduardo Fradkin, Field Theories of Condensed Matter Systems, Addison & Wesley, 
1991. 

F. Wilczek, Phys. Rev. Lett 49, 957 (1982), and Phys. Rev. Lett 48, 1144 (1982). 

D. Eliezer and G. W. Senienoff, Phys. Lett B286, 118 (1992), and references therein; 
this paper corrects the commutation relations of reference [jl9| which are inconsistent. 



G. S. Canright, S. M. Girvin, and A. Brass, Phys. Rev. Lett 63, 2295 (1989). 
Ana Lopez and Eduardo Fradkin, Phys. Rev. B47, 7080 (1993). 
D. Hofstadter, Phys. Rev. B14, 2239 (1976). 

D. A. Huse and V. Elser, Phys. Rev. Lett 60, 2531 (1988). 

A detailed derivation of this equivalence can be found in V. J. Emery, in Highly Con- 
ducting One-Dimensional Solids, edited by J. Devreese, R. Evrard and V. Van Doren 
(Plenum, New York, 1979). 

E. H. Lieb and D. C. Mattis, Mathematical Physics in One Dimension, ( Academic 
Press, New York, 1966). 

X. G. Wen, F. Wilczek and A. Zee, Phys. Rev. B39, 11413 (1989). 



[44] In this form, this theory is known as the one-component Gross- Neveu model, D. J. Gross 



and A. Neveu, Phys. Rev. DIO, 3235 (1974). We are following here reference |34 



[45] The last term in eq. (|5.12| ) is just equal to ^B'^ and it is not a full Maxwell Lagrangian 
which, instead, should look like -^{^ — B'^) where vq is the speed of propagation 
of the gauge fields. Such electric-like terms are absent at the present order but do 
get generated by renormalization effects. In any event, both the coupling constant (or 
effective "charge" ) e and the speed vq are dimensionful parameters. A full Maxwell-type 



51 



term is still irrelevant at long distances. Thus, in principle, such terms can (and will) 
be ignored. 

[46] We use the statistical mechanical sign in the definition of the /^-function. 

[47] M. P. M. den Nijs, Phys. Rev. B23, 6111 (1981). 

[48] F. D. M. Haldane, Phys. Rev. B25, 4925 (1982). 

[49] Wei Chen and Miao Li, Phys. Rev. Lett.70, 884 (1993). 

The critical value g* corresponds to a A* ~ "\/3/4 < 1. The fact that it is not precisely 
equal to one (the nominal value of the anisotropic point) is due to the approximations 
that went in taking the continuum limit. The same happens in the bosonization ap- 
proaches (based in the Wigner- Jordan transformation) in one space dimension [^. In 
any case the symmetries and, consequently, the universality class, are the same. 



52 



FIGURES 
FIG. 1. Fermion mass gaps. The full curves show the mass of the fermion species rrii = Mi/ A 

(i = 1,2) against the (normalized) coupling constant t = g/gc- The broken curves are the fermion 

masses m = M/A, predicted by the AFA. The phase transition occurs at t* = v3 where the mass 

rrii vanishes. The range of the Ising and XY phases is also shown. 
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